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We present structural properties of two-dimensional polymers as far as they can be described by 
percolation theory. The percolation threshold, critical exponents and fractal dimensions of clusters 
are determined by computer simulation and compared to the results of percolation theory. We also 
describe the dependence of the typical cluster structures on the reaction rate. 
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I. INTRODUCTION 

Polymers are known to develop various forms, for ex- 
ample linear chains, stars, rings, combs, ladders, three- 
dimensional networks and branched polymers. All these 
kinds of polymers have linear chains between the branch- 
ing point. The kind of two-dimensional polymers we 
are studying are compact sheets of monomers, which are 
linked together periodically with respect to two dimen- 
sions. Between the branching points there is no chain 
and the polymers are highly organized sheets. 

The successful "in bulk" synthesis of such two- 
dimensional polymers was reported recently j^J^ . Stupp 
et al. polymerized self-organized bilayers by two differ- 
ent so-called stitching reactions, which act within three 
distinct levels of the bilayer. Each oligomer can not have 
more than two bonds within a layer. In addition, not 
all possible bonds are present, because the stitching re- 
activity is about 30 to 50% at the upper and lower layer 
(90% within the middle layer). Is it possible to form a 
two-dimensional polymer instead of an ensemble of lad- 
der polymers? What is the minimal reaction rate, needed 
for large two-dimensional polymers? What kind of struc- 
tures will be developed at different reaction rates? 

These are questions, which can be answered by perco- 
lation theory ||] . 

II. RESULTS FOR THE STITCHING REACTION 

One of the most important results of percolation the- 
ory is related to "universality" , if the system has 
short ranged interactions. All presently available evi- 
dence strongly suggests, that the critical exponents de- 
pend on the dimensionality of the lattice only — but not 
on the lattice structure, boundary conditions and so on. 
The exponents for bond or site percolation, for square, 
triangular or honeycomb lattices etc. are the same. We 
therefore would expect, that the clusters produced with 
the above cited stitching reaction also may be described 
by these exponents. 

Does there exist a large ( "infinite" ) cluster of bonded 
sites? To answer this question, we define the fraction of 
sites in the largest cluster Poo and the probability Ps, 
that there is a "spanning cluster" of bonds, which con- 
nects two opposite boundarys of the underlying lattice. 
Then for an infinite lattice. Poo is expected to vanish for 
(P - Pc) 0+ as 

Pooip) 'X (p - pcf , (1) 

where p is the probability, that a bond is present, and 
Pc is the percolation threshold, above which an infinite 
cluster exists. With clusters of size s (number of poly- 
mers) the mean size of the finite clusters is related to the 
percolation susceptibility x (omnitting the largest ( "infi- 
nite") cluster) 



oo ' 

' (2) 

s=l 

which diverges at Pc as 

X{p) oc \p - Pcl"^ . (3) 

Finally, at p = pc the number Hs of clusters with size 
s is expected to decay with a power of s: 

ns{pc) oc s^'^ (4) 

Here /3, 7andT are the mentioned exponents. 

Stupp et al. [|| suggest, that they synthesized large 
two-dimensional polymers. Obviously, the effective reac- 
tion rate of the stitching reaction has to be larger than 
the percolation threshold pc, if one wants to synthesize 
large sheets instead of ladder polymers. Since Pc is not 
universal, we use a simple but sufficient model of the two- 
dimensional polymer, to compute Pc- In addition we will 
measure the exponents. 

As a model of the two-dimensional polymer we used a 
stack of three layers, each triangular, hexagonally shaped 
with free boundaries, as shown in figure^. At the middle 
layer, the reaction rate was 100%, i.e. all allowed bonds 
(two can react) were occupied. For that purpose, the 
bonds of the triangular lattice were choosen successively 
at random. If the number of bonds of each polymer was 
less than 2, the randomly choosen bond was occupied. 
Then we stitched the upper and lower layer using the 
same algorithm, until n+i (n_i) bonds were occupied, 
n+i {n-i) is the number of occupied bonds of the middle 
layer no times the stitching probability (reaction rate) 
p of the upper and lower layer. Note, that riQ is not 
determined by the linear lattice size P, but may vary 
slightly for each configuration. 

Because of the restriction, that each polymer may be 
stitched to at most two neighbours, we do not expect an 
infinite cluster or percolation at p = 0, i.e. without bonds 
at the upper and lower layer, but an ensemble of linear 
polymer chains. On the contrary, at p = 1 all allowed 
bonds at all three layers are occupied, resulting in an 
effective coordination number larger than two. It is worth 
noting, that even at p = 1 the probability, that there is 
only one cluster containing all polymers, vanishes with 
increasing system size L. Within this model, it is very 
unlikely to prepare very large perfect two-dimensional 
polymers. 

If we want to answer the question, if a reaction rate of 
about 40% is high enough to prepare infinite clusters, we 
have to determine the percolation threshold. As noted, 
the value of Pc is not universal, but depends strongly 
on the used lattice and bond restrictions. An unbiased 
way to estimate Pc is to plot the probability Ps of the 
occurence of a spanning cluster as a function of p for 
different lattice sizes [|. The finite size scaling relation 
for Ps simply reads 



FIG. 1. Schematic representation of the two-dimensional polymer model of linear size L = 11 and p = 0.4. 
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where P5 is the scahng function and the correlation 

length. As a consequence, different curves PPip) for 
different choices of L should intersect at p = Pc in a 
common intersection point Ps(0), apart from corrections 
to scaling and a bias by the choice of the critical expo- 
nents. Figure H shows this plot for the stitching reaction 
described above. 

From this figure we get one of our main result, the 
percolation threshold Pc = 0.110(3). Hence it follows, 
that the reaction rate of 30% to 50% reported by Stupp 
et al. |l| was high enough (within the limits of our per- 
colation model) to prepare an "infinite" two-dimensional 
polymer. 

Does the percolation of the stitching reaction lie in the 
same universality class as the standard bond percolation? 
Figure || shows a finite size scaling plot of Poo using the 
assumed exact values for (3 = 5/36 and i/ = 4/3 and 
Pc obtained from the spanning probability. Within the 
accuracy of the data, the finite size scaling assumption 
seems to be fulfilled. 

We may also compare the cluster size distribution to 
bond percolation. Using 7 = 43/18, we plotted the scal- 
ing function of the susceptibility x in figure ^. Apart from 
small deviations of the smallest system far from pc, the 
system sizes simulated seem to be in the scaling regime. 
The maximal values of x vary by orders of magnitude in 
this range of L\ 

Finally, we are interested in the geometric properties 
of the cluster structure, particularly at the percolation 
threshold Pc- At pc the largest (spanning) cluster should 
form a fractal with Hausdorff dimension dn ■ This quan- 
tity can be defined with the number of polymers N within 
a sphere of radius R and the relation 

N = const. X R'^" . (6) 

Using the data of 1000 configurations of a L = 1001 
hexagon at pc, which is shown in figure |5[ we measured 
dH{Pc) = 1.91(5) in good agreement with the result of 
percolation theory 91/48 w 1.896. 

Apart from the largest cluster, there are many smaller 
ones. Because of the statistical nature of their structure, 
we restrict ourselves to the number Us of clusters of size s 
(which is the number of polymers), which is expected to 
decrease as a power of s. Within the accuracy of the data, 
the power law seems to be fulfilled in figure Although 
the Fisher epxonent r w 1.95(10) is somewhat smaller 
than the expected one 187/91. Most probably this is 
due to the finiteness of the system. High-quality data 
[0 showed, that ris is less than the expected power law 
for small s and higher for the largest s. Therefore, we 
underestimated r systematically. 
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FIG. 2. Unbiased estimate of the percolation threshold pc = 0.110(3) using the spanning probability Ps- 
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FIG. 3. Finite size scaling of the probability P!h \ which was defined as the fraction of polymers in the largest cluster, with 
exponents v = 4/3, (3 = 5/36 and pc = 0.110(3). 
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FIG. 4. Finite size scaling of the susceptibility x'^^ with exponents v = 4/3, 7 = 43/18 and pc = 0.110(3). 




FIG. 5. Number of polymers N within a radius R of the largest cluster. The data o/lOOO configurations of L = 1001 hexagons 
were accumulated. The dashed line has the exact Hausdorff dimension (slope) dn = 91/48. 




FIG. 6. Cluster distribution of 1000 configurations of L — 1001 hexagon (751501 polymers) at the percolation threshold 
Pc « 0.110. From the slope one computes the Fisher exponents t « 1.95(10). 
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